Theory of Optical Transmission through Elhptical Nanohole Arrays 
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We present a theory which explains (in the quasistatic limit) the experimentally observed [R. 
Gordon, et al, Phys. Rev. Lett. 92, 037401 (2004)] squared dependence of the depolarization ratio 
on the aspect ratio of the holes, as well as other features of extraordinary light transition. We 
calculated the effective dielectric tensor of a metal film penetrated by elliptical cylindrical holes and 
found the extraordinarily light transmission at special frequencies related to the surface plasmon 
resonances of the composite film. We also propose to use the magnetic field for getting a strong 
polarization effect, which depends on the ratio of the cyclotron to plasmon frequencies. 

PACS numbers: 64.60.Ak, 73.23.-b; 72.80.Tm, 78.66.Sq, 77.84.Lf 
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A pioneering article of Ebbesen et. a'} reported on 
an extraordinary optical transmission through periodic 
holes array in metallic films. This was explained by the 
coupling of light with surface plasmons. In most of the 
articles^ published after Ref.[ll, the surface plasmons were 
treated as the coupled waves, propagating along both 
film surfaces (in framework of the theory described in 
Ref.i). In this approach the periodicity and spacing be- 
tween the holes were taken into account, but the shape 
of the holes was not. In Ref. 4 it was even written that 
in "the long- wavelength limit ... the transmission coef- 
ficient ... does not appreciably depend on hole-shape". 
In contrast to this, in our paper we treat the plasmon 
as the excitations localized around the holes^. It was 
shown that the present of the magnetic field transforms 
the initially circular holes into elliptical ones in the re- 
scaled virtual coordinate space. That is, it was predicted 
that the elliptical shape of inclusions should induce some 
anisotropy into the system. This prediction was appar- 
ently not known to the authors of Ref. who actually 
confirmed it in their recent experiment on light trans- 
mission through a periodical array of real (not virtual) 
elliptical holes. 

In this paper we present a theory, which in the qua- 
sistatic (long wavelength) limit explains the squared de- 
pendence of the depolarization ratio on the aspect ratio 
of the holes (experimentally observed in Ref. i7|), as well 
as other optical features. This theory is a further devel- 
opment of our methods described in Refs. dHH. We also 
note that similar effects (e.g., light polarization) can be 
reached by applying a static magnetic field. 



Let us consider a geometry which corresponds to 
the above-mentioned experiments^: a metal film with a 
square array of identical perpendicular elliptical holes. 
A monochromatic light beam of angular frequency uj im- 
pinges upon this film along the perpendicular axis y, with 
linear polarization along the principal axis x of the array 
(see Fig.[T]). 

Following Refs. [Idll, we can treat the holes as dielec- 
tric inclusions embedded in a conducting host^ii^ In this 




Z,£|| 






h 

























(b) 



FIG. 1: (a) A schematic drawing of a metal film with a 
periodic array of elliptical holes. The incident light beam 
is normal to the film surface, i.e., the ac electric field E is 
parallel to the film plane, while the wave vector is normal to 
it (i.e., k II y). (b) An elliptical hole in a unit cell when the 
main semi-axes a and c are inclined in respect to the lattice 
axes X and z by an angle qq. 



approach, the local electric potential 0("^(r) is then the 
solution of a boundary value problem based upon the 
Laplace partial differential equation 



V • EM ■ V( 



,(") 



V • OiSe ■ V0' 



(1) 



and the boundary condition (j)^"'' — Va- Here is the a- 
component of r, £/ and £m are the electrical permittivity 
tensors of the inclusions and the metal host respectively, 
Se = Em — £i, ^/(r) is the characteristic function describ- 
ing the location and the shape of the inclusions {9i = 1 
inside the inclusions and Oi = outside of them)<^ i^°d^ 

The host with the anisotropic permittivity tensor em 
can be transformed to an isotropic i'j^j using the rescaling 
of the Cartesian coordinates (^i = x/^/Exxi S,2 = 
y/ ^Eyy, ^3 = z/ y/Ezz)- Thcu cUiptic cylinder inclusion 
will be transformed into some new elliptic cylinder in the 
rescaled virtual <^-space, and the electric field E/ V0/ 
inside this single inclusion {inclined by some angle in 
respect to the main axes [see Fig. [Ijb)]} can be found 
from the system of linear equationsi^jiS. 
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FIG. 2: (Color online) The depolarization ratio of the trans- 
mitted light as a function of the aspect ratio a/c of the holes 
for different A. For the values of A from ~ 900 nm to ~ 700 
nm the curves are close enough to the square law y = (left 
dashed curve), while already for A = 550 nm and A = 500 nm 
(shown by dotted lines) the polarization ratio is essentially 
different from this law. 



0.4 




0.2 - A 



where n^p are the Cartesian components of the depolar- 
ization factor— li^. The latter can be transformed to the 
diagonal form riap = UaSap by a simple coordinate rota- 
tion r' = R{olq) -r, where ao is the angle on inclination of 
the elliptical holes from the lattice axes [see Fig. [Hb)], 
i?(ao) is the rotation matrix 

R(ao) = {(cosQfo, — sinao)i (sinao, cosao)}, (3) 

which directs the new coordinate axes along the principal 
ellipse axes. Note that naf3 (as well as Ua) now depends 
on the precise shape of the transformed inclusion, and 
therefore is a function of Sm- If the coordinate axes are 
the principal axes of the inclusion {uap = UaSap), then 
the elliptic cylindric hole of the semi-axes a and c (with 
symmetry axis along y) transforms in ^-space into a new 
elliptic cylinder with semi-axes a' — a/ ^Sxx and c' — 
cj ^Ezz for which: 



(4) 
(5) 

For simplicity, we assume the host permittivity tensor 
has the form e = eo-I+i^a, where the conductivity ten- 
sor (T is taken in the free-electron Drude approximation 
with Bq II z: 
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So is the scalar dielectric constant of the background 
ionic lattice, and / is the unit tensor. The magnetic 
field enters only through the Hall-to-Ohmic resistivity 
ratio H = pn/p = axy/cFxx = m|Bo| = ^cT, where 
LUc — eB/mc is the cyclotron frequency, r is the con- 

1/2 

ductivity relaxation time, ujp = (47re^iVo/m) is the 



FIG. 3: (Color online) Transmission spectrum through an 
elliptical nanohole arrays for the p- and s- linear light polar- 
izations. The p polarization is parallel to the [0,1] direction, 

while s is parallel to the [1,0] direction, (a) Im e'xl (Eo || 2;, 
p polarization) and Im ei^J (Eo || z, s-polarization) vs. tu/uip. 
(b) Re e'^xx and Re ei^J vs. co/cop. In both polarizations there 



is an SP resonance peak at lo 



and uj 



(shown by 



the arrow), respectively. When Eo || z, the resonance shifts 
for larger values {ujspz > i^spx) and its amplitude reduces 
drastically, (c) Transmission coefflcient T vs. frequency uj, 
for different polarizations E || 2; and E || z. The analytical 
results obtained using a dilute and MG approximations are 
indicated by dashed curves. The (connected by lines) filled 
and open symbols denote the results obtained numerically^^ 
for a square array of holes in the shape of elliptical cylin- 
ders. Inset to Fig. (a): a unit cell with the elliptic hole and 
coordinate axes. 



plasma frequency, Nq is the charge carrier concentration, 
TO is its effective mass, and /i is the Hall mobility. 

When Exx = £zz [H — 0), we have Hx — c/{a + c), 
Uz = a/(a-|-c). However, when H > the depolarization 
factor has a complicated dependence on uj. 

Solving the system of linear equations we obtain 
the following expression for the uniform electric field in- 
side the original elliptical cylinder (when Bq || z): 



E, 



7-Eo, 



(7) 



where 7 is a 3 x 3 matrix whose nonzero components are: 

"^xx — ^xx/ i^xx ^x^^xx) J ^xy — '^x^^xy / {^xx '^x^^xx^ 
IVV = 1' 7zz = £zz/{£zz - nzSSzz)- 

From Eq. ([7]) it follows that the depolarization ratio 
[i.e., the ratio of the light intensity Ix (polarized parallel 
to X-axis) to the light intensity Iz (polarized parallel to 
2;-axis), and therefore proportional to \Ex/Ez\^], can be 
written (in the case of nap = UaSa/s) as 




(8) 
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Let us consider the case of zero magnetic field, H = 0, 
when both the host and the inclusions are isotropic and, 
therefore, are characterized by the scalar tensors Em — 
£m ■ I t Si — Ei ■ I (where / is a unit matrix). We assume 
also that _Eox ~ Eqz- When cem ^ aEj and qem ^ cei, 
what might be true in the simple limit eq <C cf/lo and 
wr 3> 1 with H = (situation of Ref. I?!), then from Eq. 
(EI) it follows that 



(9) 



as in Fig. 4 of Ref. Q. 

In Fig. [2] we show the ratio ^ for different wave- 
lengths, A. Using the data^ of the complex permittivity 
Em of the evaporated gold, we can see that in the wave- 
length's range between A ^ 900 nm {em ~ —28.0 4- il.8) 
and A ^ 700 nm {em ^ —14.7 -I- il.O), the curves are 
close enough to the law {c/af' , while already for A = 500 
nm (for which em ~ —3.4 -|- zO.7), the polarization ratio 
\Ex/Ez\'^ (shown by the dotted line) is essentially differ- 
ent from this law. 

If ^ and a = c, then Ix/h — \ExxlEyy\ (see Eq. 
(HI). In Drude approximation © (and in the limit Eq 0, 
LOpT —f 1) this takes the simplest form 



Ix/Iz — \Ex/Ez 



1 - {uc/ujY 



(10) 



from which it follows that the depolarization ratio Ix/Iz 
can be made arbitrarily large, since the value lJc/^ = 
(ujcl ujp){uj I ujp)^^ = {H /ojpT){u} /uip)~^ [where H = uJcT, 
see comments to Eq. (jS])], can be made as close to 1 as 
necessary. 

Next, we approximately compute the tensor £e(^)j 
which is defined by the relation ie ■ (E(r)) — (e(r) •E(r)), 
where (...) denotes a volume average, and £(r) is the local 
dielectric tensor. In the case of the dilute collection of 
the elliptic cylinders, the tensor Eg takes the for m^'^°'^^ 
£e — — P^E • 7) where p is the volume fraction of in- 
clusions. For the case 11 = (and when Eq and the 
coordinate axes are directed along the symmetry axes of 
the ellipse), ie takes the form 



Em [1 - pSEii/{EM - niSEu)] 



(11) 



where rii is given by Eq. ([5]), and i = x, y, z. 

The frequency Wsp,j of the surface plasmon (SP) po- 
larized in the z"* direction is the one in which E [see 
Eq. ([7]) and Eq. (fTT|) ] becomes very large even for a 
very small applied field. This condition is satisfied when 
^M;iiiysp,i) — n'iSEii{uJsp,i) = 0. Substituting Eq. © into 



this, and letting ujpT 



one obtains 



(12) 



For Ei ~ 1 this simplifies into uJsp,x = Wp-^/l — Ux- 

In Maxwell-Garnett (MG) or Clausius-Mossotti 
approximatioiiii, expressions for Ef. and Ugp can be ob- 
tained directly from Eqs. pTjl . p^ . obtained in dilute ap- 
proximation just by formal substitution rii ni(l— p)^ 
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FIG. 4; (Color online) The same as Fig. |3jc), but vs. wave- 
length A. The peak at A ~ 600 — 700 nm corresjjonds to 
Eq W X (i.e., p-) polarization (shown in Fig. 2 of Ref. while 
the peak at A ~ 400 nm corresponds to iJo || polarization 
(not shown in Fig. 2 of Ref. [7j). Up = 2, ■ 10^^ rad/s. 



When the aspect ratio c/a tends to infinity the con- 
sidered geometry transforms into the case of the parallel 
slabs, which can be solved exactly: l/ce = pm/em + 
Pi/ej. Using the form ([6]), we can find the resonance fre- 
quency LOsp = oJp\Jpi / {pmEi +Pi)- This coincides with 
Eq. (|12p in the limit Ux ^ 1- The other exact solvable 
geometry, parallel cylinders, for which 0, does not 

give any resonance: e^ = Pm^m+Pi^i — 1 — PA/(^p/^)^- 

When the cylindrical holes are arranged on a two- 
dimensional periodic lattice^'^ a more suitable approach 
is a Fourier expansion technique^". Since Oj{y) and 
■)/)(") = — r*^"-' are now periodical functions, they 
can be expanded in a Fourier series. This transforms 
Eq. ^ into an infinite set of linear algebraic equations 

for the Fourier coefficients V'g;"'' — ^ /y V'*-"' (i")^"**''"'^!^, 
where g = {2'K/d)(mx,my,mz) is a vector of the appro- 
priate reciprocal lattice, rrii are the arbitrary integers, d 
is a lattice constant, and V is the volume of a unit cell. 
After solving a truncated, finite subset of those equa- 
tions, we can use those Fourier coefficients V'g"'' o^'" 
der to calculate the bulk effective macroscopic electric 
permittivity tensor Eg, using the procedure described in 
Ref. 5,10-11,. Note that the Fourier coefficient 6g of the 
(?/(r)-function of the elliptical hole (inclined by the angle 
ao in respect to the main lattice axes) has the form^iiS 
= {ATrac/d^g±) Ji(gj_)[sin(|5;,|/i/2)/|g^|], where Ji(x) 
is a Bessel function and g±{a)= [a^{gx cos a+ Qz sina)^-!- 
c^{gz cos a— (/a; sina)^]^/^. 

In Figs.[3]Ja) and [3Kb) we show the imaginary and real 
parts of Ee{oj) vs. lojiOp, respectively. The curves with- 

(e) 

out points are the two principal in-plain components Exx 
and E^zz (corresponding to polarizations Eo||a; and EqU^, 
respectively) as obtained in the dilute and MG approx- 
imation. The full and open squares in Fig. [DJa) denote 
the same quantities, but now for a square lattice (of lat- 
tice constant d) of elliptical holes with the same aspect 
ratio and volume fraction p as in MG and dilute approx- 
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imations. In this case, Sxx and efz are calculated by the 
Fourier expansion technique^ mentioned above. 

Finally, in Fig. EJc), wc show the calculated trans- 
mission coefficient T for the dielectric functions shown 
in Figs. [3l^a) and [3l^b). The dependence of the trans- 
mission coefficient T on frequency w, [as well as on 
the wave length A (see Fig. |4|)] for different film thick- 
nesses, can be obtained from the effective value Se using 
the known expressiorti^ for T = where d = {1 — 
rl^)/[e M-ix )-rl2(iMiX% and na = {l-N)/{l + N), 
N = ^ee(tj), X — [uj/cjhN — {lo / uip){u)ph/ c)N , and h is 
the film thicknessi^ The transmission coefficient T(ijj) [see 
Fig. IHl^c)] shows the characteristic "extraordinary trans- 
mission" peaks expected, based on Figs. |31[a),(b). Since 
the angular dependence of the permittivity tensor in the 
rotated coordinate system e' = R{a + aQ)-e-R^^{a + aQ) 
[where R{a -I- ao) is given by Eq. ([3]) and a is the po- 
larization angle i.e., the angle between the vector Eq 
and the axes of the ellipse, see Fig. [TJb)] is a cosine-like 

[e.g., {e^xly = e-xl cos^{a -I- ao) -f elz' sin^(Q; + ao)], the 
angular dependence of the transmission coefficient T(a) 
looks (depending on the values of a/c and uh/c) also as 
a cosine-like. This explains the Malus' lawi^ of the light 
transmission, T, observed in Refs. 

The amplitudes and frequencies of the peaks (for both 
polarizations Eq || x and Eq || z) depend on the aspect- 
ration a/c. For polarization along z axis the resonance 
frequency ujsp z shifts for larger value and its amplitude 
reduces drastically. Estimating the imaginary part of 
fe-^^ = e-^^ — £a/ [see Eq. (fTTj) ] at the resonance fre- 
quency uJspi, we found that the ratio Im de^xli^spx) /^'^ 



5eiz{ujspz) is of the order ~ (riaj/n^)^/^ = (c/a)^/^. For 
the aspect ratio a/c = 0.3 this is of the order ~ 20, in 
agreement with our numerical calculations [see Fig.[3]Ja)]. 
One might expect that in the extraordinary light tran- 
sition the maximum at tUspz will also be much smaller 
when compared to the maximum at cosp x , but [as we can 
see in Figs. [SJc) andH], they are of the same order. In 
Fig. |4]we show the data, presented in Fig. ^c), in terms 
of wavelength A. The calculations were performed for the 
systems with aspect ratio a/c = 0.3, volume fraction of 
holes, p = 0.031, ej — 1, ujpT ~ 40, ujp ~ 3 ■ 10^^ rad/s 
(a typical for Au valued, so that ujph/c = 1), iVpT = 40, 
^ = ujph/c = 1, and the film thickness h — lOOnm. 

In summary, we have studied analytically and numer- 
ically the extraordinary transmission through perforated 
metal films with elliptical holes. We have explained an- 
alytically the optical features found experimentally and 
described in Ref. 0- Our numerical results are in good 
agreement with experimental data. We also propose to 
use the magnetic field for getting a strong polarization 
effect, which depends on the ratio lo^/^p- As a material 
which may be suitable for this purpose the bismuth can 
be considered, where the low carrier density (~ 3 x 10^^ 
cm^^) can make the carrier cyclotron energies, cjc, equal 
to or greater than the plasmon energy, ujp)^ 
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